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ABSTRACT 

We present stellar metallicities in Leo I, Leo II, IG 1613, and Phoenix dwarf galaxies derived from 
medium (F390M) and broad (F555W, F8I4W) band photometry using the Wide Field Camera 3 
(WFC3) instrument aboard the Hubble Space Telescope. We measured metallicity distribution func¬ 
tions (MDFs) in two ways, I) matching stars to isochrones in color-color diagrams, and 2) solving for 
the best linear combination of synthetic populations to match the observed color-color diagram. The 
synthetic technique reduces the effect of photometric scatter, and produces MDFs 30-50% narrower 
than the MDFs produced from individually matched stars. We fit the synthetic and individual MDFs 
to analytical chemical evolution models (CEM) to quantify the enrichment and the effect of gas flows 
within the galaxies. Additionally, we measure stellar metallicity gradients in Leo I and 11. For IG 1613 
and Phoenix our data do not have the radial extent to confirm a metallicity gradient for either galaxy. 

We find the MDF of Leo I (dwarf spheroidal) to be very peaked with a steep metal rich cutoff 
and an extended metal poor tail, while Leo II (dwarf spheroidal), Phoenix (dwarf transition) and IG 
1613 (dwarf irregular) have wider, less peaked MDFs than Leo 1. A simple CEM is not the best fit 
for any of our galaxies, therefore we also fit the ‘Best Accretion Model’ of Lynden-Bell 1975. For 
Leo II, IC 1613 and Phoenix we find similar accretion parameters for the CEM, even though they all 
have different effective yields, masses, star formation histories and morphologies. We suggest that the 
dynamical history of a galaxy is reflected in the MDF, where broad MDFs are seen in galaxies that 
have chemically evolved in relative isolation and narrowly peaked MDFs are seen in galaxies that have 
experienced more complicated dynamical interactions concurrent with their chemical evolution. 
Subject headings: galaxies: abundances; galaxies: dwarf; galaxies: evolution; (galaxies:) Local Group 


1. INTRODUCTION 

Dwarf galaxies are important constituents of the uni¬ 
verse because they are both the most numerous type of 
galaxy and test beds for examining galactic evolution 
on small scales. Additionally, dwarf galaxies show vary¬ 
ing and extended perio ds of star formation (SF ) over 
the age of the universe (|Weisz et al.ll201lH 1201411 . One 
method of examining these populations is through metal¬ 
licities, specifically through metallicity distribution func¬ 
tions (MDF). Metallicity measurements of large numbers 
of stars over the full range of metallicities are crucial in 
studying the populations of dwarf galaxies because they 
allow us to construct robust MDFs, examine any struc¬ 
ture indicating sub-components, and measure gradients 
across the galaxy. The star formation history (SFH), 
accretion, outflows via supernova and stellar winds, and 
galaxy interactions that cause tidal or ram-pressure strip¬ 
ping all play a part in shaping the MDF. 

The shape of the MDF offers clues to the galaxy’s 
evolution, which can be characterized by fitting them 
to analytic chemical evolution models. The most ba¬ 
sic analytic chemical evolution model is a simple closed 
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box, where gas turns into stars, the stars evolve and 
return enriched gas to the interstellar medium (ISM), 
which is again formed into stars. In this model no ma¬ 
terial enters or leaves the system. However, the simple 
model is often an inaccurate description of the chemi¬ 
cal evolution in dwarf galaxies because the simple model 
overpredicts the number of metal poor stars (i.e. the 
G dwar f problem). Semi-analytic chemical evolution 
mode ls (iLanfranchi &: Matteuccil I20I0I: lHendricks_et_^ 
2Q14) and hydrodynamic sim ulations (|Marcolini et al. 
2006 . [20?)M iReyaz et al.ll20n9[ l utilize the galactic SFHs 
of dwarf galaxies, and nucleosynthetic yields from type I 
and H supernoya (SN), to predict the gas inflow and out¬ 
flow necessary to match abundance patterns from mul¬ 
tiple elements and the oyerall MDFs. In lieu of these 
complex models the simple analytic chemical eyolution 
models offer a concise way to quantify and compare the 
MDFs of different galaxies, proyiding insight into the en¬ 
richment history of galaxies, especially when only mod¬ 
eling the oyerall MDF of the dwarf galaxy without infor¬ 
mation on a abundances. 


Recent works have examined the stellar metallicities 
found in dwarf galaxies by using medium and high reso- 


lution spectra on 8-10 meter telescopes (iTolstov et al. 

20031 IShetrone et all I2001L 120031 

20091 Kirby et al. 

20IIbl 20131 Starkenburff et al. 12013 

: Koch et al.l 120131 

Lemasle et al.ll2014HHendricks et al.||2O140. These works 


examine variations in a elements, abundance patterns in 
the lowest metallicity stars, abundance patterns in r and 
s process elements, and the overall MDFs to better un¬ 
derstand the processes that are important in the chemical 
evolution; processes that include the SFHs, the IMF, stel- 
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TABLE 1 

Observational data 


Galaxy 

Name 

Observation 

dates 

Total Exposure Time (sec) 
F390M F555W F814W 

Leo I 

2011, Jun 02 & Mar 25 

21,024 

1,760 

1,500 

Leo H 

2012, Mar 24 & 25 

10,080 

1,760 

1,500 


2013, Mar 30 

10,080 

— 

— 

IC 1613 

2011, Dec 16 k. 20 

15,720 

1,200 

1,224 

Phoenix 

2012, Jan 30 & 31 

16,320 

1,300 

1,340 


lar nucleosynthesis, and the time-scales for the formation 
of galaxies. One of the largest samples of stellar metallici- 
ties (> 3000 stars) comes from iKirbv et al.l ()2011bl[201311 
who present spectroscopic MDFs for 15 dSph and 7 dirs. 
They also fit their MDFs to analytical chemical evolu¬ 
tion models, extended the stellar-mass stellar-metallicity 
relation and determined that the MDF shapes vary de¬ 
pending on the morphology such that dSph tend to have 
narrower MDFs than dIrs. 

While individual spectra of stars in dwarf galaxies pro¬ 
vide abundance and kinematic information, they are still 
difficult to obtain. Spectra of stars in more distant Local 
Group (LG) objects are limited to small sample sizes by 
the long exposure times and large telescopes required to 
make the observations. While spectral targets are lim¬ 
ited to the few brightest stars, photometric metallicities 
probe deeper and thus sample stars that are typically 
more representative of the metallicity distribution. Pho¬ 
tometric metallicities, though not as accurate as spectra, 
provide measurements for every star in the field, allowing 
us to increase samples sizes by an order of a magnitude. 
This allows us to probe galaxies farther out in the LG, 
especially when using the resolving power of the Hub¬ 
ble Space Telescope (HST). We obtained HST images 
in metallicity sensitive filters (F390M), and 2 wide band 
filters (F555W and F814W) for four LG galaxies. By 
measuring every star in the field we can build up larger 
samples of stellar metallicities in these dwarf galaxies. 

In this work we present photometric metallicities of 
individual stars in four LG dwarf galaxies: Leo I, Leo 
II, IG 1613 and Phoenix. We chose galaxies that span 
different morphologies, masses, SFHs and distances from 
the Milky Way (MW). Section 2 describes our two meth¬ 
ods of measuring the photometric MDFs, including a new 
synthetic color-color diagram method, which is similar to 
the synthetic color magnitude diagram (GMD) method of 
deriving SFHs. In section 3 we compare our metallicities 
to literature values, includi ng recent work with over lap- 
ping spectroscopic targets (jKirbv et al.ll2nilbll2013jl . In 
Section 4 we fit the MDFs with chemical evolution mod¬ 
els to quantify the enrichment of the galaxy. In section 5 
we present metallicity gradients from the central regions 
of each galaxy. In section 6 we discuss implications of 
our results on theories of galaxy evolution, and conclude 
in section 7. 

2. OBSERVATION 

The WFC3 observations were obtained between 
March, 2011 and Marc h 2013 as part of HST proposal 
12304 ([Holtzmanl 1200811 . The dates and exposure times 
of the observations are given in Table |TJ We used the 
reduced images from the STScI pipeline, which performs 
bias, flat-field, and image distortion corrections. We ad¬ 


ditionally used a charge transfer efficiency (GTE) cor¬ 
rection module provided by STScI 0. Magnitudes were 
measured using the point spread function fitting photom¬ 
etry packa ge DOLPHOT , which is a modified version of 
HSTphot ijPolphinl 1200011 . The photometric catalog was 
culled to reject objects based on goodness of fit and pro¬ 
file sharpness as recommended by the DOLPHOT man- 
uaE 

We adopt reddeni ng and distance values reported in 
iMcGonnachiel (|20I2ll to calculate the absolute magni¬ 
tudes. All magnitudes are reported in the Vegamag sys¬ 
tem. Table [2] lists some basic observable quantities for 
these dwarf galaxies. 

3. METHODS 
3.1. Deriving Metallicity 

The general technique to measure photometric metal¬ 
licities relates color to metallicity. The stellar proper¬ 
ties that control observed colors are metallicity, effective 
temperature and surface gravity. However, effective tem¬ 
perature and surface gravity are dictated by the mass, 
metallicity and the current evolutionary stage of a given 
star. For populations of comparable age the color is di¬ 
rectly related to the metallicity. However, an individ¬ 
ual giant branch star can be redder either because it has 
higher metallicity or because it is older, leading to a color 
degeneracy between age and metallicity. For a mixed 
aged population, the younger giants are hotter and more 
massive than older giants of the same metallicity so the 
color-metallicity relation no longer holds. 

Using specifically designed HST WFC3 filters which 
separate the effects of metallicity and temperature on 
color, we can break the age-metallicity degeneracy. 
The metallicity sensitivity of the (F390M-F555W) color 
comes from the F390M medium band filter which cov¬ 
ers the Ga H & K spectral features, one of the strongest 
metal absorption features in the visible spectrum. The 
temperature sensitive color, (F555W-F8I4W), uses two 
broad band hlters which cover mostly continuum. 

The color-color diagram (F390M-F555W, F555W- 
F8I4W), has been shown to effectively separate the com- 
peting color chan ges due to metallicity and temperature 
(iRoss et aril20I4ll . breaking the age-metallicity degener¬ 
acy. In the color-color diagram the color change between 
a 12 and 4 Gyr isochrone of the same metallicity is on 
order of a few hundredths of a magnitude, (demonstrated 
in Figure 5 of Ross et al. 2014). Therefore, these two 
colors can be used to measure individual stellar metallic¬ 
ities of populations of mixed ages and abundances, like 
dwarf galaxies. 

3.2. Measuring Individual Stellar Metallicities 

We used only giant branch stars with errors less than 
0.05 mag in all filters. These error cuts resulted in sam¬ 
ples of 3,449 stars, 444 stars, 896 and 578 stars for Leo 
I, Leo H, IG 1613 and Phoenix respectively. 

Stellar metallicities are assigned in the (F390M- 
F555W, F555W-F8I4W) color-color diagram where each 

® We used the alpha version (2013) for the CTE em¬ 
pirical pixel based corrections for WFC3/UVIS CTE located 
http: //www.stsci.edu/hst/wfc3/ins_performance/CTE/ 

° The 2014 DOLPHOT 2.0 WFC3 module manual can be found 
at http://americano.dolphinsim.com/dolphot/ 
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TABLE 2 

Dwarf galaxy properties 


Galaxy 

Ra 

Dec 

D 

(kpc) 


rii 

(’). (pc) 

M* 

(IO^Mq) 

^Total 

[WMq) 

(T* 

(km s“^) 

Mhi 

( 1 O®M 0 ) 

17 HI 

(km s“^) 

type 

IC 1613 

01 04 47.8 s 

+02 07 04 

721 ± 5 

-14.6 

6.81, 1496 

100 

79.5 


65 

25 

dir 

Phoenix 

01 51 06.3 

-44 26 41 

406 ± 13 

- 10.1 

3.75, 454 

0.77 

3.3 


0.12 

10 

dTrans 

Leo I 

10 08 28.1 s 

+12 18 23 

254 ± 17 

-11.9 

3.40, 251 

5.5 

2.2 

9.2 

N/A 

N/A 

dSph 

Leo II 

11 13 28.8 

+22 09 06 

233 ± 15 

-9.8 

2.60, 176 

0.74 

0.97 

6.6 

N/A 

N/A 

dSph 


Note. — Data are taken from pTolstov et al.l II2009I) and IMcConnachiel II2012IL IMcConnachiel (120121) reported the stellar mass, M*, assuming a 
stellar mass-to-light ratio of 1. Virial masses reported from [Mated il998D . 


lO 

LO 

lO 


o 

CD 

n 


1.0 1.2 1.4 i.e 

F555W-F814W 


1.8 1.0 1.2 1.4 1.6 1.8 1.0 1.2 1.4 1.6 1.8 1.0 1.2 1.4 1.6 1.8 

F555W-F814W F555W-F814W F555W-F814W 




[M/H] [M/H] [M/H] 


Fig. 1.— The top panels show the color-color diagrams of giant branch stars in M92, NGC 6791, NGC 104 and NGC 5927. The middle 
panels show Hess diagrams of the residuals from the observations minus the best combination of synthetic populations normalized to the 
square root of the synthetic populations. The bottom panels show the MDFs from the synthetic method in solid black as well as the MDF 
from individually matched stars, as dashed lines. 


giant bran ch star is matched to the closest Dartmouth 
isochrone (|Dotter et al.ll200§l in a grid (spaced by 0.05 
dex in [Fe/H] and assuming solar [ct/Fe]). The isochrones 
were empirically corrected to align with globular clus¬ 
ters of known metallicities (M92, NGC 6752, NGC 104, 
NGC 5927 and NGC 6791 with [Fe/H] = -2.30,-1.45, 
—0.70, —0.40 and +0.40 respectively). Th e empirical cor¬ 
rections were derived in iRoss et al.l (|2014fl using the same 
HST hlters used in this study. This method of deriv¬ 
ing photometric stellar metallicities from HST photome¬ 
try was tested against the listed globular clusters and 
shown to produce metallicities with errors of 0.15-0.3 
dex , with larger erro rs occurring at lower metallicities. 
See lRoss et al.l (120141 1 for the full details on the metallic- 
ity derivations from colors. 

Although each giant branch star is matched to an 
isochrone grid of [Fe/H] assuming a solar abundance ra¬ 
tio, many stars are known to have a enhancement, espe¬ 
cially low metallicity ones. The colors we measure are a 
product of [Fe/H] and the a enhancement; therefore, it 
is more accurate to think of our measurement as an indi¬ 
cator of the total metallicity, [M/H], that represents the 
intrinsic combination of [Fe/H] and [a/Fe], even though 


we are matching each star to an isochrone grid of [Fe/H]. 
Total metallicity, [M/H], and the iron abundance, [Fe/H], 
are often used interchangeably, because they are usu¬ 
ally very similar; only in the cases of large abundance 
variatio ns (e.g a e n hance ments) do the two values dif¬ 
fer. In iRoss et all (120141) we quantified the amount of 
color change expected from a enhancements and provide 
a means of calculating the [Fe/H] if the a enhancement 
is known. 

3.3. Synthetic Color-Color Diagram Method of deriving 

MDFs 

In addition to measuring metallicities of individual 
stars based on color, we have also developed a synthetic 
color-color diagram technique to derive a MDF from pho¬ 
tometry. The technique uses isochrones and an initial 
mass function (IMF) to make a model color-color dia¬ 
gram of the synthetic populations in order to reproduce 
the observed color-color diagram. 

The technique is similar to the synthetic GMD method 
used to derive SFHs (see Tolstoy et al. 2009). In the 
GMD the position and number density of stars depend 
upon the IMF, age, and metallicity. As pointed out by 
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iDolphinl (|2002f) . the CMD of any complex stellar system 
will be a composite of all the individual stellar popu¬ 
lations that comprise the system. Therefore it is pos¬ 
sible to model a large number of simple stellar popula¬ 
tions (SSPs) in various combinations to reproduce the 
observed data, including the observational errors. The 
most likely combination of populations of various ages 
and metallicities that best match the observations will 
give the SFH. 

In the color-color diagram the position and number 
density of stars depend upon metallicity and the IMF of 
a population, and is independent of age, unlike a CMD. 
The synthetic color-color diagram method of deriving 
MDFs solves for the weights of a given linear combina¬ 
tion of populations that equals the observed color-color 
diagram; the weights from the linear combination give 
the metallicity distribution of the population. The ad¬ 
vantage is that a given population will occupy a locus 
where color and number density of stars found at each 
location within the color-colo r diagram is dictat ed by 
stellar evolution and the IMF (|Tolstov et al.ll20d^ . 

To create the synt hetic color-color diagrams a Kroupa 
IMF ([Kroupal l2001f ) and colors from the Dartmouth 
isochrones are used to initially populate a grid of Hess 
diagrams of SSPs. The isochrone colors were empir- 
ically corr e cted f ollowing the calibration described in 
iRoss et ^ (|2014fl . The grid spans a metallicity range 
—2.5 < [Fe/H] < -1-0.5 in steps of 0.05 dex. Each SSP 
Hess diagram was blurred using the measured photomet¬ 
ric errors from each dwarf galaxy. The linear combina¬ 
tion of SSPs that reproduces the observed color-color di¬ 
agram are calculated from the array of synthetic Hess 
diagrams. 

The isochrone color spacing and the bin spacing of the 
Hess diagrams (0.04 by 0.04 mag) require that the metal¬ 
licity spacing (and subsequent color change) be larger 
than the Hess bins. In some regions of the color-color di¬ 
agram isochrones spaced 0.05 dex of [Fe/H] apart have a 
color change less than 0.04 mag. If the color change is less 
than the bin size only one of the weights for the linear 
combination of isochrones will be positive the rest will 
have negative and nonphysical weights. Additionally, 
in the synthetic color-color diagram, when the synthetic 
SSPs are blurred by the photometric errors (up to ^0.07 
mag near the bottom of the giant branch), isochrones 
closer than 0.15 dex produce nonphysical negative solu¬ 
tions. To account for the issues arising from the isochrone 
color spacing and the bin spacing in the Hess diagrams, 
we combine the weights from various metallicity combi¬ 
nations to ensure we sample all metallicities equally. We 
start with the full range of metallicities (-2.5 < [Fe/H] < 
0.5, and spacing of 0.05 dex), calculate the linear com¬ 
bination, eliminate the metallicities that have negative, 
nonphysical solutions and recalculate the linear weights. 
This leads to uneven spacing in metallicity due to the 
smaller color separations at low metallicity and larger 
color spacing at high metallicity. We repeat the calu- 
lation with different metallicity spacings (O.I, 0.15, 0.2, 
and 0.25 dex) and shifts of 0.05 dex in order to fill in the 
metallicity spacings. The weights from various metal¬ 
licity combinations are combine to ensure we sample all 
metallicities equally. 

The MDF weights are checked by creating a synthetic 
color-color diagram to compare to the observations. We 


examine the residuals between our observed Hess dia¬ 
gram and a synthetic Hess diagram to ensure that > 99% 
of the residuals deviate by less than Scr of the overall 
residual within the Hess diagram. The residual Hess di¬ 
agram is created using the following equation: 

. , , ,, observed Hess — synthetic Hess 

residual Hess diagram = - , - 

V synthetic Hess 

( 1 ) 

To test the synthetic color-color diagram method of 
deriving MDFs we performed this analysis on glob¬ 
ular clusters of known metallicity, specifically M92, 
NGC 6791, NGC 104, and NGG 5927, with [Fe/H] = 
—2.3,—1.45,—0.70,—0.40 respectively. For these GCs 
we ad opted all of the parameters reported in IRoss et ^ 
(|2014f) . where the individual MDFs were derived. Here 
after we will refer to the MDFs derived using the syn¬ 
thetic color-color diagram method as ‘synthetic MDFs’, 
and the MDFs created from individually measured 
metallicities as ‘individual MDFs’. The synthetic MDFs 
recovered peak metallicities that are all within 0.05 dex 
of the literature values. Additionally we find the syn¬ 
thetic distributions to be over two times as narrow as 
the distributions found from fitting each star individu¬ 
ally (cr =0.15, 0.19, 0.24, and 0.19 for M92, NGG 6752, 
NGC 104 and NGC 5927, respectively). The narrower 
MDFs are expected becanse this method accounts for 
the photometric errors, while the individually measured 
metallicities do not. 

Figure [T] shows the process and final results of the 
synthetic color-color diagram method of deriving MDFs. 
The top panels display the observed seqnence of giant 
branch stars in the color-color diagram. The middle pan¬ 
els show the Hess diagram of the residuals as described 
in Equation (1). The bottom panels show the resulting 
MDFs from the synthetic method (solid line) as well as 
the individual MDFs (dashed line). The MDFs in the 
bottom panels illustrate the utility of using the synthetic 
method over the individual star method, wherein the syn¬ 
thetic MDFs are 30 - 50 % narrower than the individual 
MDFs. The synthetic method is statistical in nature, 
therefore it does not provide individual metallicities for 
each star, but rather gives the relative number of stars 
at each metallicity for the population as a whole. 

4. COMPARING METALLICITIES 
4.1. Star by star metallicity comparison 

We directly compare the subset of our individually 
measured stellar sa mple that ov e rlap with spe ctroscopic 
me asurements f r om iKirbv et al.l (l2011bL l2013f) in Figure 
m IKirbv et al.l (l2011bL 1201311 used medium resolution 
spectra and spectral synthesis of Fe I absorption lines to 
measure metallicities for 814, 256, and 125 stars in Leo 
1, Leo 11, and 1C 1613, respectively. Due to the greater 
spatial extent of their sample we were only able to com¬ 
pare a total of 170 stars; 108 from Leo 1, 54 from Leo 11, 
and 8 from 1C 1613. The standard deviation of metallic¬ 
ity differences is 0.38 dex, a wider spread (b y 0.1 to 0.2 
dex) t han found in the globular clusters from IRoss et al.l 
(|2014f) . The larger spread between photometric and spec¬ 
troscopic metallicity can be partially attributed to the 
difference in photometric errors between the dwarf galax¬ 
ies and the globnlar clusters. The photometric errors 
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along the dwarf galaxies giant branch are 2 to 4 times 
larger than equivalent absolute magnitudes found in the 
globular clusters, which translates to a (5[Fe/H] of ~0.1 
dex. 

The larger spread might al so be partia l ly acc ounted for 
by a abundance variations. iRoss et all (j2014ll find that 
variations in a abundance cause color changes analogous 
to metallicity changes of a few tenths of a dex (specifically 
A[Fe/H]/Aa ^ 0.65 — 0.34 across the metallicity range). 
Without knowing the intrinsic a abundance, and assum¬ 
ing it to be solar, naturally our metallicity measurement 
will be lower than the actual [Fe/H]. In general dwarf 
galaxies tend to show less a en hancements than glo bular 
clusters of similar metallicity. iKirbv et al.l (|2011ar) find 
that the a abundance distributions slowly evolve from 
large a enhancement at low metallicity to roughly solar 
ratios at [Fe/H] close to a tenth solar. Any a enhance¬ 
ment will cause an underestimate of the metallicity using 
the photometric method. 

In addition to the greater difference between photo¬ 
metric and spectroscopic metallicities, the photometric 
metallicities deviate more towards lower metallicity. This 
is not unexpected as we assign metallicities with an 
isochrone grid assuming no a enhancements. For any star 
that has an enhanced a abundance the assigned metallic¬ 
ity from isochrones will be lower than the intrinsic stellar 
metallicity because the color of a star of solar a abun¬ 
dance will have the same color as an a enhanced star 
with lower [Fe/H]. 

The stars with photometric metallicities measured to 
be [M/H] < —2 tend to be the most discrepant, with 
differences > 1 dex. The main difference is due to the 
decreasing color change as a function of decreasing metal¬ 
licity, which means that a small random variation in color 
for a bluer (metal poor) star will produce a larger discrep¬ 
ancy in the reported metallicity than the same amount 
of variation in a redder (metal rich) star. 



Spectroscopic [Fe/H] 

Fig. 2.— The comparison of photometric and spectroscopic 
metallicities for 108 stars from Leo I, 54 from Leo II, and 8 from 
IC 1613 shows that the photometric metallicities deviate to the 
metal poor side and the deviations are worse at lower metallicities . 
Spectroscopic metallicities measured bv IKirbv et al.l (I2011bl . l2013ll . 


4.2. MDFs comparisons 

We used the synthetic color-color diagram technique 
to derive MDFs for the four dwarf galaxies in addition 
to the individually measured metallicities that were com¬ 
piled to make MDFs. Figure [31 following the same layout 
as Figure in shows the results from the synthetic MDF 
derivation method. The top panels show the color-color 
diagrams for the four dwarfs, the middle panels show 
the residual Hess diagrams, and the bottom panels show 
the resulting synthetic (solid line) and individual (dash- 
dotted line) MDFs. The individual and synthetic MDFs 
both show similar shapes and peaks, although the syn¬ 
thetic MDFs are narrower, as was noted in section 3.3 
from the MDFs of globular clusters. 

For the dwarf galaxies the average metallicities of the 
synthetic MDFs are O.I to 0.2 dex more metal rich than 
the individual MDFs, although the MDF peaks from the 
two methods are more closely aligned. Additionally, the 
globular cluster MDFs from the two methods are not 
offset. This leads us to believe that the underlying pop¬ 
ulation distribution is causing the offset. For the dwarf 
galaxies the difference in average metallicity from the two 
methods can be attributed to the differences in the metal 
poor tails of the MDFs, where the individual MDF has 
photometric errors that propigate into a wider spread in 
metallicity. Additionally, the fact that the color change 
due to metallicity is smaller at low metallicities means 
that for the same amount of color error there will be a 
larger spread in metallicities at low values. The large 
extent of the metal poor tail in the individual MDF is 
the main driver in the offset of the average metallicities 
for the 2 methods. 

Despite the offset in average metallicity we believe 
the synthetic method provides a more accurate MDF 
shape than the MDF produced by individually measur¬ 
ing metallicities with isochrone matching. The synthetic 
method includes stellar sequence information from the 
color-color diagram, and it systematically accounts for 
the photometric errors, whereas the individual method 
relies only on the photometry. 

In Figure |3| we compare our two methods of deriv- 
ing MDFs t o the s pectr oscopic MDFs available from 
IKirbv et al.l (|20IIbl I2013D . Table [3] reports the average 
metallicities, widths, and number of stars measured us¬ 
ing the three different MDF derivation meth ods for each 
galaxy. The s pectroscopic metallicities from IKirbv et al.l 
(|20IIblI20I3I1 are more spatially extended than our pho¬ 
tometric data. Since we are probing a smaller area within 
each galaxy, it is expected that our MDFs will have some 
differences compared to the spectroscopic MDFs. The re¬ 
sults for each dwarf galaxy are discussed in the following 
paragraphs. 

For Leo I we find <[M/B]individuai > = -1.43. The 
synthetic MDF is narrow (cr = 0.16), mostly Gaussian 
with an extension on the metal poor end, and shows an 
abrupt cutoff on the metal rich side of the distribution. 
The shapes, peaks and widths are consistent with metal¬ 
licities found in the literature. IKirbv et al.l (I2011bll found 
the median [Fe/H] = —1.42 {a = 0.331. iBosler et'al] 
(l2007fl measured a MDF (for 101 RGB stars) peaked at 
[Ga/H] = —1.34 (tr = 0.21 dex), using the infrared Ga- 
triplet method of measuring st ellar metallicities. Using 
a different Ga-triplet calibration IGullieuszik et al.l (|2009f) 
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Fig. 3.— The top panels show the color-color diagrams of giant branch stars in Leo I, Leo II, IC 1613 and Phoenix. The middle panels 
show Hess diagrams of the residuals from the observations minus the best combination of synthetic populations normalized to the square 
root of the synthetic populations. The bottom panels show the MDFs from the synthetic method as the solid histogram, while the dashed 
lines are the MDFs derived from star individually matched to isochrones. 


found [Fe/H] ~ —1.37 (a = 0.18 dex) using 54 stars. 

For Leo II we find <[M/H]i„diaid«ai > = -1.77 with 
the MDF showing a more extended tail at low metallici- 
ties. The peak value and overall MDF shap e are consis¬ 
tent wi th those reported in the literature. iKirbv et al.l 
(|20IIbll measured a MDF with a peak [Fe/H] = —1.71 
from the medium resolution spectra of 256 stars. A 
near IR photome t ric st udy of RGB stars in Leo II by 
iGullieuszik et al.1 (j2008[ l found a MDF peak, [M/H] = 
-1.64, when they account for the mean age (9 Gyr). 
iBosler et al.l (120071) measured stellar metallicities of 74 
stars using the infrared Ca-triplet of RGB stars, the 
MDF they measured p e aked a t [Ca/H] = —1.65 (a = 
0.17 dex). iKoch et ^ (|2007ari also measured spectro¬ 
scopic metallicities using the Ga H infrared method and 
found the mean metallicity of Leo H to be [Fe/H] = —1.73 
from 52 stars. 

For IG 1613 we find <\M/Yi]individuai > = -0.99, 
which is more metal rich than the average s pectroscopic 
metallicity measured bv IKirbv et al.l (|2013[) . <[Fe/H]> 
= —1.19. Anot h er ave rage metallicity measurement from 
iSkillman et al.l (|2014D finds <[M/H]> = —1.3 from a 
SFH study of a non-central field. The discrepancy be¬ 
tween the higher average metallicit y that we measure 
com pared to the va l ues r eported bv IKirbv et al.l (|2013[) 
and ISkillman et aP (|2014[) . is partially accounted for by 
the fact that we are looking at different regions of the 
galaxy, and thus a younger population of stars. SFHs of 
IG 1613 show the star forming regions mo ve continuously 
inwar d over the lifetime of the galaxy (iBernard et al.l 
I2007D . where SF is associated with higher average metal¬ 
licity. Additionally, ISkillman et al.l (I2003D find the stellar 
metallicities of IG 1613 to progress from [Fe/H] = —1.3 
to —0.7 over the age of the galaxy. Our field of view 
looks at the inner ^5 % of the galaxy, while the spec¬ 
troscopic metallicities are from stars spanning the whole 


galaxy and the SFH comes from a field 5.5’ from the 
center. While the average metallicity we find is more 
metal rich than most studies, we are more metal poor 
than the [Fe/H] = —0.67 from the high resolution stud y 
based upon 3 supergiants by (iTautvaisiene et all 120071) . 
We believe the average metallicity we measure is reason¬ 
able for the young central region of IG 1613 that our field 
of view covers 5%) even though it is more metal poor 
than the average galactic metallicity measured in other 
studies. 



Fig. 4.— The synthetic MDFs (black lines) are compared to the 
individual MDFs (blue lines). For three of the galaxies spec tro- 
scopic MDFs (red lines) as measu red bv|Ki rbv et ^ J20liljl are 
also shown as a comparison. The IKirbv et alTT20^^) MDFs are 
true [Fe/H], as opposed to our measurements of [M/H], which are 
a combination of [Fe/H] and [o/Fe]. 
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For Phoenix we find <[M/Y{\individuai > = -1.17. 
Upon first examination this average metallicity is in 
conflict with metallicities reported in the literature, al¬ 
though none of the literature values are spectroscopi¬ 
cally derived. Most of the re ported metallicit i es are; from 
SFH studies, like the one bv iHoltzman et aP (I2000D who 
found <[Fe/H ] > ^ —1.7. The radial SFH study by 
iHidalgo et al.l (|2013ll found the old population to have 
Fe/H] ~ —1.7, and the younger populations to have 
Fe/H ~ -1 .2. Additional studies find [Fe/H] = —1.9 , 
—1.81±0.10 (jOrtolani &: Grattonlll988HHelcl et al.lll999H 
A study of variable RR Lyrae stars in Phoenix predicts 
<[Fe/H]> = —1.55 or —1 .75 based on two dif ferent pe¬ 
riod metallicity relations (|Ordonez et al.ll2014fl. Accord - 
ing to the mass-metallicity relationship (|Lee et al.ll2006H . 
Phoenix, which is the same stellar mass as Leo 11, should 
have a similar metallicity ([Fe/H] ^ —1.7); yet, the av¬ 
erage metallicity we measure is higher than the galactic 
average metallicity reported in other studies. 

However, the field of view for our observations covers 
the innermost region 15%) of Phoenix, where SF has 
occurred as recent as 100 Myr ago and higher metallic¬ 
ities are expected due to the continuous star formation 
occurring in t he central region. Ad ditionally, the radial 
SFH study bv IHidalgo et al.l (I2013D find the central re¬ 
gions of Phoenix to have a mean metallicity of ~ —1.2, 
which is in agreement with our results. As with IC 1613, 
we believe the higher average metallicity we measure re¬ 
flects the greater amount of enrichment that occurs in 
the central regions of dwarf galaxies. 


5. ANALYTIC MODELS OF CHEMICAL EVOLUTION 

The chemical enrichment within galaxies provides in¬ 
formation on the gas accretion, gas expulsion, the inter¬ 
action history and star formation history. Matching an¬ 
alytic models of chemical evolution to observations will 
constrain different evolutionary scenarios. It stands to 
reason that dwarf galaxies could be represented by simple 
models due to their isolation and relatively less compli¬ 
cated dynamics. Therefore we begin by fitting a simple 
model of chemical evolution to our observed MDFs. 

The simple model is essentially a leaky box, where gas 
is allowed to leave the galaxy, but no accretion is occur¬ 
ring. Additionally we assume p is the effective yield, not 
the true nucleosynthetic yield, where the effective yield 
is a product of both the amount of metals created, and 
the amount blown out of the galaxy through winds. The 
effective yield p is in units of Zq . The shape of the MDF 
for the simple model follows the functional form as de¬ 
fined bv iPagell (I1997D : 


dN 

d[Fe/H] 



10[Fe/H] 

p 


( 2 ) 


The simple model has a tendency to overproduce metal 
poor stars, whether in dwarf galaxies or the solar neigh¬ 
borhood. We find that the simple model is often un¬ 
able to completely reproduce the M DF of dwarf galax¬ 
ies, in agreement with the findings of (|Kirbv et al.ll20lfbl 
l2013f) . It has been theorized that the paucity of metal 
poor stars can be explained by gas infall at later times 
([Prantzo'i l2008f) . Therefore we also tested an analyti¬ 
cal model which includes inf alling gas. W e fit th e ‘Best 
Accretion Model’ defined by iLvnden-Belll (|1975f) which 


allows for the gas mass to start small (or at zero), rise 
to a maximum, then approach zero again as all the mass 
is accumulated in stars. The accreted gas is assumed to 
be metal free. Lynden-Bell found a relation between the 
gas mass, g, stellar mass, s, and total final mass, M, that 
permits an analytic solution to the differential metallicity 
distribution (also see Pagel 1997), defined as: 

where all quantities are in units of the initial gas mass. 
When M = 1 the equation reduces to a closed box where 
g = 1 — s. When M is larger than 1 it essentially be¬ 
comes a measure of the total amount of gas that has 
entered the system over the lifetime of the galaxy so M 
can be thought of as an accretion parameter. However, it 
should be noted that while Equation (3) produces the de¬ 
sired MDF shapes, there is no physical rationale behind 
the actual form of the equation. Equation (3) simply 
assumes a quadratic relation between gas mass and the 
stellar mass that peaks then decays as a function of stel¬ 
lar mn^. The MDF of the accretion model, as defined 
by IPagell (1199711 , follows the form: 

dN lOPe/H] 1 -k S(1 - jA) 

d|Fe/H]“ p 

(4) 

where s must be solved for numerically from the equa¬ 
tion: 

[Fe/H](s) = (5) 



For both the simple and the accretion model, the 
metallicity peak of the model directly increases with p, 
as long as instantaneous mixing of the ISM is assumed. 
The analytic form of the differential metallicity distri¬ 
bution for these models requires the instantaneous re¬ 
cycling and mixing approximations, where all SN yields 
are immediately and uniformly returned to the ISM, re¬ 
spectively. The instantaneous recycling approximation 
does not reproduce the characteristic patterns seen in a 
element abundances. However a abundances cannot be 
determined with photometric data alone, therefore we 
could not constrain more sophisticated models that ac¬ 
count for time dependent recycling. The analytic models 
used in this work neglect some of the physics known to 
be important for galactic chemical evolution, however, 
these models provide insight into the difference in evolu¬ 
tion between galaxies. 

5.1. MDF Truncation of the Metal Rich End 

For dSph, one leading hypothesis is that their chem¬ 
ical evolution was interrupted by the removal of gas. 
This follows from the morphology-density relationship 
where most dwarf galaxies close to a large primary 
are dSph, and devoid of gas (|Grcevich fc PutmanI 120091 : 
iSoekkens et al.l[20I4f) . Environmental mechanisms such 
as tidal stirring and ram pressure stripping have been 
invoked to transform close satellites from a star forming. 




















































TABLE 3 
MDF PROPERTIES 


Galaxy Individual MDF Synthetic MDF Spectroscopic MDF 

Name <[M/H]> a Num. of stars <[M/H]> a Num. of stars <[Fe/H]> a Num. of stars 


Leo I 

-1.43 

0.43 

3,449 

-1.34 

0.21 

4,296 

-1.45 

0.32 

814 

Leo II 

-1.77 

0.72 

444 

-1.54 

0.27 

780 

-1.63 

0.40 

256 

IC 1613 

-0.99 

0.54 

896 

-0.84 

0.30 

2,111 

-1.19 

0.37 

125 

Phoenix 

-1.17 

0.67 

578 

-1.14 

0.28 

1,399 

— 

— 

— 


gas rich, dir, into a quiescent, gas p oor dSph that are 
found close to a much larger primary (iLin h FabedllQ^ 
Maver et al.l 120011: iGrebel et al.ll^03HKazantzidis et al l 


201 Ih . iMaver et alJ ( 20061) used simulations to shows 


that gas rich dwarfs can be stripped of their gas over 
a few pericentric passages as they travel through the hot 
halo of a larger primary. Ram pressure stripping effec¬ 
tively cuts off all star formation, stopping the chemical 
evolution prematurely, likely producing a steep trunca- 
ti on on the m e tal ric h end of the MDF. 

IKirbv et al.l (|2013ll show that a sharp metal rich cut¬ 
off of the MDF is primarily seen in dSph type galaxies, 
not dirs. Additionally, they added a ram pressure strip¬ 
ping term to the simple chemical evolution model to fit 
the metal rich cutoff, however, this model was found to 
not have the best match to the dSphs in their sample. 
Since the accretion models fit all the dSph galaxies bet¬ 
ter than the simple models we would like to incorporate 
a metallicity cutoff into the accretion models, that would 
simulate ram pressure stripping. 

A more basic approach to testing a metal rich cutoff 
is achieved by changing the upper metallicity bounds to 
a more metal poor value when integrating the equations 
of chemical evolution. We ran the simple and accretion 
chemical evolution models with an imposed metal rich 
cutoff, truncating the CEM at a lower metallicity than 
the model would naturally evolve to, in order to approx¬ 
imate the effects of ram pressure stripping. A sharp 
cutoff such as this is somewhat nonphysical, since ram 
pressure stripping does not immediately remove all gas, 
additionally the SF would not be immediately cutoff at 
a particular metallicity across the galaxy. However, this 
simplistic approach allows us to quickly and efficiently 
test our models for signs of truncation. 

Besides the dSph galaxies, dir are expected to have a 
steeper MDF on the metal rich side than the chemical 
evolution models predict. The MDF of a star forming 
galaxy should not look like a completed chemical evolu¬ 
tion model, rather, the MDF should show a more abrupt 
cutoff on the metal rich side because dir are still forming 
stars and enriching their ISM and by definition are not 
at the end of their chemical evolution. Fitting CEM with 
an imposed metal rich cutoff will also simulate an incom¬ 
plete chemical evolution, in addition to a CEM that has 
been halted due to ram pressure stripping. 


5.2. Best fit models 

To hnd the best model parameters we perform 10^ tri¬ 
als of the simple model while varying the effective yield, 
p. We calculate the least squares value between each 
model and MDF to determine the best parameters for 
the individual and synthetic MDFs for each galaxy. We 
also calculate the two sided confidence interval to find 
the 1 a range of model parameters. For the Accretion 


Model we preform 10'^ trials while varying the effective 
yield, p, and the extra gas parameter, M. For the ac¬ 
cretion + truncation models we run 10® trials varying 
the effective yield, p, the extra gas parameter, M, and 
the cutoff metallicity. The best fit parameters for each 
model and MDF are listed in Table S) 

To fit the individual MDFs the three chemical evolu¬ 
tion models have been convolved with a Gaussian with 
a dispersion equal to the dispersion of the difference 
between the photometric and spectroscopic metallici- 
ties (cr=0.38 dex) to account for the spreading of the 
MDF from the photometric errors. The synthetic MDF 
method models the photometric errors, therefore we do 
not convolve the chemical evolution models with any ad¬ 
ditional errors when fitting to the synthetic MDFs. 

The individual MDFs do not show an abrupt metal 
rich cutoff, mostly due to the scatter in metallicity that 
comes from the photometric technique as can be seen in 
Figure El However, in the narrower synthetic MDFs we 
do see truncation of the MDFs on the metal rich side, 
indicative of ram pressure stripping. In the accretion 
model the larger the M parameter the more the CEM 
takes on a Gaussian shape. While this is good for fit¬ 
ting the narrow peak of the synthetic MDFs, it does not 
account for the asymmetry of the extended metal poor 
tail and the sharper truncation of the MDF on the metal 
rich side as shown in Figure El The accretion -|- trunca¬ 
tion model allows for a narrower MDF peak, and keeps 
the asymmetry by truncating the metal rich side of the 
distributions, which produces a better model fit for the 
synthetic MDFs in all four dwarf galaxies, as shown in 
Figure [71 


5.2.1. Leo I 

Leo I is one of the more distant MW companion dwarf 
spheroidal galaxies (254 kpc, McConnachie 2012); addi¬ 
tionally it is receding quickly from the MW (Kad=167.9 
km s“^, Vtan =101 km s“^, Sohn et al. 2013), making 
it one of the more un usual dSph galaxies. Since its dis¬ 
covery in the 1950’s (jHarrington fc WilsonI [T950ll Leo I 
has been studied by many groups to better understand 
the SFH, dynamic properties, metallicities and chemi¬ 
cal enrichment history. Here we compare our chemical 
evolution models to results from previous studies. 

We find the accretion -|- truncation model is the best 
fit model for the synthetic MDF. The synthetic MDFs 
shows a narrow peak at [M/H] ~ —1.35, a slight metal 
poor tail, and an abrupt cutoff on the metal rich side 
of the distribution. The simple model underpredicts the 
number of stars in the peak, and overpredicts the amount 
of metal poor stars. The accretion model hts the MDF 
better than the simple model, however this model does 
not ht the truncation on the metal rich side of the MDF, 
and subsequently overpredicts the number of metal rich 
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[M/H] [M/H] 


Fig. 5.— The individual MDFs and Poisson uncertainties are 
shown with the best fit simple and accretion chemical evolution 
model. The number or stars in each MDF are listed for each galaxy. 


stars. We find that the narrowness of the synthetic MDF, 
and the metal rich truncation is best fit by the accretion 
+ truncation models, with a large accretion parameter 
of M = 6 . 60 l 2 ^g, yield of p = 0 . 071 q;qo 5 and a cutoff 
metallicity of [ Fe/H] = —1.0 . 

iKirbv etldl (j2011bL l201f{[l also modeled the MDF of 
Leo I, with a simple, pre-enriched and accretion model. 
They found the accretion model was the best fit, with 
M = 7.9 and p = 0.043, which is within errors bnt 
smaller than our accretion model parameter of M = 
and the yield we derive, p = 0.068, is higher than 
the one they derive, but our MDF has a higher average 
metallicity (see discussion i n Sec tion 4.2). 

iLanfranchi fc Matteuccil (|2010f) created detailed chem¬ 
ical evolution models for Leo I and II, that use nucleosyn¬ 
thesis for type I and II SN, an exponentially decreasing 
accretion parameter, the SFH, and a abundances and 
the MDF as inputs to solve for the outflowing wind, and 
star formation efficiency. The detailed chemical evolu¬ 
tion model that best matches the Leo I abundances pre¬ 
dicts that the galaxy started with a low star formation 
efficiency (0.6 Gyr“^) for the first 5 Gyr, then had an¬ 
other 7 Gyr episode of SF starti ng at 9 Gyr with substan¬ 


tial g alactic winds throughout ( Lanfranchi fc Matteucci 


The evolution model of ILanfranchi fc Matteucci 
agrees with the color magnitude diagram (GMD) 
based SFHs showing that Leo I had slow ancient SF, with 
an increase in the SFR around 7-8 Gyr ago, and another 
increase to its highest rate 2-3 Gyr ago, a fter which the 
galaxy stopped forming stars I Gyr ago (iGallart et al 


Caputo et al.llI99ftlHernandez et al1l200Cll:iDolphin 


1999; 


2002 ; 


Weisz et al1l20I4ll. 


Fig. 6 .— the synthetic MDFs with Poisson uncertanties and the 
best fit simple and accretion chemical evolution models are shown. 
The Poisson uncertainties are calulated from the total number of 
observed stars modeled by the synthetic color-color diagrams. 

most likely pysing into the MW halo for the first time 
(iLmine et al.l[2?)Tlf ). 

The synthetic MDF of Leo II is best fit by the accre¬ 
tion -|- truncation model with an accretion parameter of 
M = 4 . 012 ^ 4 , yield of p = 0.038tQ;go8 and a metallicity 
cutoff of [Fe/H] = —I.O. The simple model underpredicts 
the number of stars found in the MDF peak, and over¬ 
predicts the number stars in the metal poor tail. The 
accretion model fits the peak and metal poor tail, but 
over predicts the number of stars on the metal rich end 
of the distribution, explaining why the best fit model 
includes the metallicity truncation term. 

Our chem i cal evolution results are in agreement with 
IKirbv et al.l (|2011bL[20I3ll who found the accretion model 
to be the best fit for Leo II, with an accretion pa¬ 
rameter of M = 3 . 11 q '5 and yield of p = 0.028. For 
our accretion model we measured M = 4 . 1 ^ 3 ;? and 
p = 0.038. ILanfranchi fc Matteuccil (l201flf l modeled the 
detailed chemical evolution of Leo II based on the overall 
MDF and the a elem ent abundance ratios measured by 
iShetrone et al.l (|2009ll . Their best fit model predicts one 
long (7 Gyr) episode of SF starting at 14 Gyr, with a 
lower SF efficiency than seen in Leo I (0.3 Gyr“^) with 
very hig h wind efficiency throughout. T he chemical evo¬ 
lution of iLanfranchi fc Matteuc^ (|2010ll is in good agree¬ 
ment with th e SFH derived from the w ide field photomet¬ 
ric survey of iKomivama et ^ (|2007| 1. They found that 
Leo H evolved with a low SF rate up to 8 Gyr ago when 
SF stopped in the outer regions, from 8 to 4 Gyr ago the 
central star forming region continually shrank until SF 
essentially stopped. 


5 . 2 . 2 . Leo II 

Leo H is also one of the more distant MW companion 
dwarf spheroidal galaxies (233 kpc, McGonnachie 2012). 
Unlike Leo I, the proper motion study of Leo H reveals 
the galactocentric velocity is mostly tangential {vtan = 
265.2 kms~^), with only 21.5 kms~^ in the radial direc¬ 
tion. Subsequent dynamic modeling of the large veloc¬ 
ity and large distance of Leo H suggest that Leo H is 


5 . 2 . 3 . IC 1613: 

As one of the nearest gas rich dwarf irregular galax¬ 
ies at a distance of 721 kpc there have b een numerous 
studi es of the gas and stars in IG 1613 (jMcGonnachid 
l2012f) . IC 1613 is similarly distant from M31 (517 kpc) 
however it is not considered to be a satellite of either the 
MW or M31 ([McConnachiel [2M^ . but without data on 
the proper motion the membership of IC 1613 to either 
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TABLE 4 

Chemical evolution model parameters 


Galaxy 

Name 

MDF type 

Simple Model 

P 

Accretion Model 
p M 

Accretion + truncation 
p M [Fe/H] cutoff 

Best Model 

Leo I 

Leo II 

IC 1613 

Phoenix 

Indiv 

Indiv 

Indiv 

Indiv 

0.054 t°oZl 
0.038 ZZl 

n 09ft +0.041 

U.ZZ6 _o .045 

Q 1 -| Q +0.015 
U.ilo o.oi 2 

0.054 

0.038 

n 1 su +0-034 
0-Io4 _0.032 
r» 1 1 9 +0.020 
U.ilZ o.oi4 

8 .2/1;® 

3.0/];® 

l-Sias 

l-°-0.5 

0.050 

0.038 ZZl 
n 1 fto +0.036 
U.184 -0.032 

0 112 

u.Ti^ -0.018 

7.7//® 

3.0//,® 

^•^-0.5 

-0.10 

-0.40 

-0.30 

-0.20 

A-l-T 

S 

A-l-T 

A 

Leo I 

Leo II 

IC 1613 

Phoenix 

Synth 

Synth 

Synth 

Synth 

u.uo^ _0 QIO 

u.uuo _o.oi8 
0.182 

0.086 

0.068 

0.038 ZZl 
0.196 

0.092 

y.o_2 0 

^•^-2.1 

4 1+2-0 
^•^-2.1 
o 4+2.0 

o.^_2.4 

0.070 ZZl 
0.038 t°oZl 

0 200 

u.zuu -0.036 

u.uy^ -0.018 

a c+2.2 
6-612.8 

^•'^-2.4 

o '7+2.0 

' -1.9 

o q+2.6 
'^•'^-1.6 

-1.0 

-1.0 

-0.40 

-0.60 

A-l-T 

A-l-T 

A-l-T 

A-l-T 


is not definitive. Various H I stu dies have measured the 
total H I mass as 3-6 xIO^Mq (jLake fc SkillmanlllQSOl : 
iSilich et al.ll2006l l. Additionally, the studies show com¬ 
plicated morphology including numerous H I arcs, holes 
and shells produced from the ongoing SF within IC 1613 
as evidenc ed by the presence of OB associations and H 
II regions (iGarcia et al.ll2010ll . 

For IC 1613 we find the accretion -|- truncation model 
with an accretion parameter of M = 3 . 7 tu 9 : a yield of 
p = 0.2010 03 ® and a cut off metallicity of [Fe/H] = —0.40 
is the best fit to the synthetic MDF. IC 1613’s MDF 
shows a broad peak at [M/H] = —1.00, an asymmet¬ 
ric Gaussian shape with an extended metal poor side of 
the distribution. The simple model underpredicts the 
number of stars in the peak, overpredicts the metal poor 
tail, and slightly overpredicts the number of stars on the 
metal rich side of the MDF. The difference between the 
accretion model and the accretion -|- truncation model is 
small, and almost entirely due to the metal rich side of 
the MDF. Since IC 1613 is an isolated galaxy that has 
a substantial gas mass and shows no sign of interaction, 
ram pressure stripping cannot be the process truncat¬ 
ing the MDF. Instead, the steeper metal rich side of the 
MDF reflects the continuing chemical enrichment of IC 
1613. 

Comparing our work with iKirbv et al.l (|2013f) we find 
that our accretion model with M = 3.7l® g and p = 
0 . 20^0 036 has a similar accretion parameter but a higher 
effective yield compared to their M = d.dtui ^-nd yield 
of p = 0.075. The higher effective yield we measure 
is potentially due to the different regions of the galaxy 
sampled. We cover a small 5%) ce ntral field while 
the metallicities from IKirbv et al.l (|2013[ 1 cover the entire 
galaxy. Eve n though our accr etion models find similar 
parameters, IKirbv et al.l (j2013l l find that a pre-enriched 
model has a more likely fit than the accretion model. 
The pre-enriched model is a generalization of a closed 
box model that starts with some initial metallicity and 
does not accrete additional gas over its lifetime. 

5.2.4. Phoenix 

Phoenix is a transition type galaxy at a distance of 415 
kpc a nd an absolute V magnitude of -9.9 (jMcConnachiel 
UHl). There is no H I emission detected within Phoenix 
itself, however, H I gas has been detected at 4 .5’ to 9’ 
from the center of the galaxy (jYoung et al.ll200'i1l indicat¬ 
ing recent gas expulsion. Young et al. also state that the 



[M/H] 


Fig. 7.— The synthetic MDFs with the Poisson uncertainties 
are shown with the best fit simple + truncation and accretion + 
truncation chemical evolution models. The Poisson uncertainties 
are calulated from the total number of observed stars modeled by 
the synthetic color-color diagrams. 

expelled gas is likely linked with the most recent episode 
of SF (~100 Myr ago; Bianchi et al. 2012). 

We find the accretion -|- truncation model to be the 
best ht model for the MDF, with an accretion param¬ 
eter of M = 3.3tf 0 , a yield of p = 0.094+g Q 3 g, and a 
cut off metallicity of [Fe/H] = —0.6. The simple model 
underpredicts the number of stars in the MDF peak and 
overpredicts the amount of metal poor stars in the tail 
of the distribution. The difference between the accretion 
model and the accretion -|- truncation is small; both mod¬ 
els slightly overpredict the number of metal poor stars 
while on the metal rich side of the MDF the accretion + 
truncation model has a slightly better agreement. 

Ground-based pho to metry of Phoenix by 
iMartmez-Delgado et al.l (1199911 noted the presence 
of a structurally distinct inner population with mostly 
young stars in the eastwest direction, while isophots 
of the older component are aligned in the north-south 
direction. Another ground based wide-held (26 x 26 
arcmin^) photometric study shows the young (< 1 
Gyr) inner population to have a disk like distribution 
([Battaglia et al.l[2012bll . While the accretion -|- trunca¬ 
tion model is found to be the best ht, evidence for ram 
pressure stripping is not seen in any of the population 
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studies, nor would it be expected in an isolated galaxy. 

I n a study of the radi al SFH properties of dwarf galax¬ 
ies iHidalgo et al.l (|2013ll additionally modeled the chem¬ 
ical evolution for Phoenix. Their best fit model predicts 
an initial episode of SF lasting 3-4 Gyr with very lit¬ 
tle or no metal enrichment that accounts for 75% of 
the total cumulative stellar mass. Additionally, the SF 
episode is longest in the central regions with shorter SF 
episodes in the outer regions. To account for the mini¬ 
mal amount of enrichment ([Fe/H] = —1.67) in this first 
phase their model predicts a high SFR that would ex¬ 
pel most of the metals formed in this episode before 
they can mix with the interstellar medium. The second 
phase of SF produces most of the metal enrichment, and 
mainly occurs in the central regions where the gas den¬ 
sity is still sufficiently high to support continued SF, they 
find the metallicity of this burst to be [Fe/H] = —1.08. 
We find the MDF peak to be [M/H] = —1.17, which 
is con sistent with the enrichment values IHidalgo et al.l 
(|2013ll found in the central regions of Phoenix from 
their radia l SFH study. Howev er, much like the initial 
metallicity IHidalgo et al.l (|2013ll found, other photomet¬ 
ric studies measur ed the overall metallic ity in Phoenix as 
<[Fe/H]> ^ -1.7, (iHoltzman et al.ll2000H . [Fe/H] = —1.9, 
( Ortolani fc GrattonI I1988D . <[Fe/H]>= —1.81 ± 0.10 
( Heldet ayi999ll . 

6 . METALLICITY GRADIENTS 

Metallicity gradients found in stellar populations offer 
clues about past SF, as well as interaction history since 
mergers are thought to remove gradients. Many dwarf 
galaxies have detectable metallicity gradients (e.g. For¬ 
nax, Getus, Garina, Sculptor, Sextans, Tucana, Leo I, 
Leo H, and Draco, Andromeda I-HI, V and VI. From 
Harbeck et al. 2001; Tolstoy et al. 2004; Battaglia et al 
2006; 2011, 2012; Monelli et al 2012; Kirby et al. 2011), 
while others like Ursa Minor and Ganes Venatici do not 
show measurable gradients (Kirby et al. 2011). 

Metallicity gradients can be produced by increased SF 
(and thus enrichment) in the central regions due to in¬ 
creased gas supply there. The gas density increases with 
the gravitational potential and thus accumulates towards 
the centers of galaxies. Additionally, star formation is 
proportional to the gas density thus producing increased 
SF and enrichment towards the insides of the galaxy. 

Environmental interactions are expected to modify the 
stellar distributions, diluting any gradients , and possi¬ 
bly in ducing morphological transformations (iLokas et al.l 
l2012f) . Additionally, recent simulations of the chemo- 
dynamical evolution of dwarf galaxies have shown tidal 
interactions to be a n efficient mechani sm to remove 
metallicity gradients (iNichols et al.l 1201411 . Also, obser¬ 
vations by IHidalgo et al.l (|201,3[1 ti^e shown that radial 
SFHs and metallicity gradients are consistent with the 
SF in dwarf galaxies being quenched in the outer regions 
due to limited gas availability. 

Gradients in age and metallicity have been seen for a 
wide range of physical characteristics of dwarf galaxies. 
The trend of younger and higher metallicity stars increas¬ 
ing towards the centers of galaxies spans many galactic 
properties such as total mass, gas mass, metallicity, ve¬ 
locity dispersion and environment. The fact that gra¬ 
dients are consistently found even across a wide range 
of physical characteristics hints that population gradi¬ 


TABLE 5 

Metallicity gradients 


Galaxy 

name 

d[M/R]/d6 
(dex deg~^) 

d[M/R]/dr 
(dex kpc”^) 

d[M/B]/d{r/rh) 
(dex per ) 

Leo I 

-1.52 

-0.34 

-0.086 ±0.04 

Leo II 

-11.55 

-3.02 

-0.54 ±0.10 

IC 1613 

-2.48 

-0.20 

-0.29 ± 0.29 

Phoenix 

-1.84 

-0.26 

-0.11 ± 0.20 


ents are intrinsic to dwarf galaxy formation and evolu¬ 
tion. However, at large radii, many of these differences 
tend to disappear, suggesting some physical processes 
may not affect all galactocentric radii equally (e.g. the 
UV background could affect the inner and outer regions 
differently). 


Leol 



Leoll 

d[M/H]/d(r/r,^)= -0.53 



IC1 613 

dlM/Hl/dfr/rJ- -0.29 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 


^ / % 

Fig. 8. — The gradients are in units of dex of [M/H] per r;,,, where 
r^ is the 1/2 light radius as defined by IMcConnachid 11201211 . 

The radial metallicity gradients we measured are re- 
ported in Tabled! and shown in FigurelH We find metal¬ 
licity gradients in Leo I and Leo H. In Phoenix and IG 
1613 the errors on the gradients are consistent with no 
slope. Much of the error is due to the small radial ex¬ 
tent the WFC3 field of view covers in these galaxies, thus 
we cannot confirm or rule out a metallicity gradient in 
Phoenix or IC 1613. 

In Leo I age and metallicity gradients have been seen 
in various star formation histor ies, variable star p opula- 
tions, and metallicity studies. iHeld et ahl (|2010f) found 
evidence for a radial gradient in th e 1-3 Gyr stellar 
popul ation using near IR photometry. iGullieuszik et al] 
(l2009fl found a shallow metallicity gradient of -0.27 dex 
kpc“^ as part of their metallicity study of Leo I, which is 
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shallower than the gradient we meas ure, -0.34 de x kpc~^ , 
in the same units. In the work by iKirbv et al.l (l2011br) 
gradients for Leo I are reported as (i[Fe/H]/fi(r/rc) = 
—0.11 dex per Tc, while we find a similar gradient for 
Leo I of —0.09 ± 0.07 dex per r^. Where the core radius 
(rc=240 pc) or the half light radius (r/i=251 pc) both 
give the same gradient value to 2 significant figures. Our 
ph otometry covers Leo I out to ^ 0.7 rh, while the work 
by IKirbv et al.l (l2011bf) extends ^4 times farther for Leo 
1 . 

In Leo IL the wid e field photometric survey by 
iKomivama et al.l (120071) found evidence of a radial gradi¬ 
ent in the HB morphology, where red HB stars are more 
ce ntrally concentrated than blue HB stars. In the work 
bv iKirbv et al.l (|2011bD the gradient for Leo II is reported 
as (i[Fe/H]/d(r/rc) —0.21 dex per Tc, while we find a 
steeper gradient of —0.54 ± 0.10 dex per r/j. Where the 
core radius (rc=180 pc) and the 1/2 light radius (r?i=176 
pc) both give the same gradient value to 2 significant fig¬ 
ures. Our photom etry covers Leo II o ut to about 0.6 rh, 
while the work bv IKirbv et al.l (l2011bD extends ^3 times 
farther. 

In IC 1613 the wide field photometric study of 
iBernard et al.l (j2007D looked at radial differences in the 
CMDs and concluded that there is an age gradient, how¬ 
ever, they could not confirm a metallicity gradient in IC 
1613 from their broadband ground based photometry. 

The radia l SFH study of Phoenix bv iHidalgo et al.l 
(l2?)n9l IWl measured the radial SF in Phoenix while 
simultaneously deriving metallicities. They found that 
Phoenix had longer periods of SF in the central regions, 
and decreasing average metallicities as function of ra¬ 
dius, with the central regions having an average [Fe/H] 
^ —1.4, decreasing outwards to —1.7, although they do 
not quantify a metallicity gradient for the galaxy. 

7. DISCUSSION 

7.1. Morphological transitions and implications from the 
MDFs and chemical evolution models 

As the name suggests, a dwarf Transition galaxy is 
considered to be in a state of morp hological transition 
between a dir and dSph type galaxy (iLin fc Fabeilll983l: 
Maver et al.ll200ll : iGrebel et ai.ll2003HKazantzidis et al l 
20111) . This idea is supported by the presence of the 
morphology-density relationship, where dSph are mainly 
found in dense environments close to a much larger 
primary while dIr and dTrans are found in isolation 
(|Weisz et al.ll2011bl) . The dwarf galaxy morphology seg¬ 
regation has led to an environment dependent explana¬ 
tion invoking tidal stirring and ram pressure stripping to 
transform a dIr close to a larger galaxy into a dTrans 
then a dSph over subsequent orbital passages around 
the primary (iMayer et al.l 120071 : iPeharrubia et al.l 120081 : 
iKazantzidis et al.ll201lD . 

However, some dwarf galaxy properties conflict with 
the morphological transition scenario. DTrans have sim- 
i lar gas fraction s (jWeisz et al.l l2011aD , and UV fluxes 
(jLee et al.ll2011l ) as dirs while the Ha emission in dTrans 
is markedly less and the typical stellar mass of dTrans 
is ^ 3-4 times lower than that of dIrs. DTrans are 
mostly located in isolated regions (like dirs), but dTrans 
luminosities , like d Sphs, are lower than dirs. This led 
IWeisz et al.l (|2011aD to suggest that dTrans appear to be 


lower mass versions of dirs that lack recent SF, and not 
an evolutionary bridge between dSphs and dirs. Addi¬ 
tionally, should dTrans be progenitors of dSph then they 
should be seen in similar locations as both dSph and dirs, 
but they are found only in the same regions as dirs. 

Further evidence against a morpho logical transition 
theory comes from IKirbv et al.l (I2013D who found sys¬ 
tematic differences in the MDF shapes of dirs and dSphs. 
They found the dSph MDF shapes were systematically 
narrower and more peaked while dirs tend to have wider 
MDFs that more closely resemble a simple chemical evo¬ 
lution model. The narrower dSph MDF is best fit by 
chemical evolution models with large amounts of accre¬ 
tion which is at odds with the proposed mechanism to 
transform dirs into dSphs, i.e. gas stripping. 

Accretion has been shown to be an important compo- 
nent in dwarf galaxy evolution (jBrook et al.l 20141). Also, 
the a nalytical ‘Best Accretion’ model of lLvnden-Bel]l 
(IT^ imply that dSph tend to require large accre- 
tion parameters (M values) to fit the observed MDFs 
(IKirbv et al.l 120131) . However, the morphology-density 
relation points to gas stripping and tidal disruption as 
important factors in shaping dSph. Simulations have also 
shown that gas stripping and tidal disrupti on is integral 
in cre ating many of the the dSph seen today (jMaver et al.l 
I2006D. Further e xamin ation of the Best Accretion model 
of iLvnden-Belll (I1975D shows that there is no physical 
basis behind the form of Equation (3) other than that it 
produces the desired MDF shapes. Equation (3) simply 
assumes a non-linear relation between gas mass and the 
stellar mass that follows a quadratic relation that peaks 
then decays as a function of stellar mass. Due to the con¬ 
tradictory implications from the accretion model and the 
gas stripping caused by gravitational interactions known 
to be integral in shaping dSphs, in addition to the lack of 
physical meaning behind the functional form of the ac¬ 
cretion model we suggest that the interpretation of the 
accretion portion of the accretion model should be left 
for more sophisticated chemical evolution models. In lieu 
of drawing interpretations from the models themselves, 
we can use the differences in the model parameters in the 
context of the dynamic histories of each galaxy to inform 
our discussion. 

Examining the dSph type galaxies in our sample (Leo 
I and H), we find that the MDF of Leo H is more similar 
to the MDFs from IC 1613 (dir) and Phoenix (dTrans) 
than to Leo I’s. Leo H is also one of the more distant 
dSph associated with the MW (233 kpc), although it i s 
slightly closer than Leo I 125 4 kpc) (lMcConnachiell2012D . 
Dynamical studies of Leo I I (iKoch et al.ll2007bD and wide 
field photometric studies (|Coleman et al.ll2007D show no 
significant signs of tidal distortion. iLeoine et al.l (j201lD 
studied Leo IPs proper motion and found nearly all of 
its velocity is in a tangential component {vtan = 265.2 
km s“^) with only a small radial component (vrad = 
21.5 km s“^). ILeoine et al.l (|2011l) take this to mean 
that Leo H either has a low-eccentricity orbit, or is near 
perigalacticon or apo galacticon, and the l ack of evidence 
for ti dal disruptions (jKoch et al.l l2007bl : iColeman et al.l 
I2007D likely rules out the latter. 

Currently the only known mechanism to transform a 
rotating dir into a dSph is through gravitatio nal influ¬ 
ence from a larger primary (iMaver et al.ll200H) . Conse- 
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quently, it can be assumed that Leo II has experienced 
some tidal interactions with the MW in the past. How¬ 
ever, Leo II shows no signs of dynamical interactions, 
beside the fact that it is a dSph. Moreover the similar¬ 
ity of the MDF of Leo II, IC 1613 and Phoenix suggest 
these galaxies experienced comparable conditions during 
their chemical evolution. Since IC 1613 and Phoenix are 
isolated and evolved without outside influence, it could 
be inferred that Leo II also chemically evolved without 
dynamical influence from the MW. The SFH of Leo II 
shows that 95% of the stars were formed over 6.3 Gyr ago 
(jWeisz et al.ll20I-ll : with so many of the stars formed at 
early times Leo II has had ample time to be transformed 
from a disk to a spheroid after most of the chemical evo¬ 
lution occurred. We do find the accretion -|- truncation 
model to be the best fit for Leo II, where there the cutoff 
on the metal-rich side of Leo IPs MDF is shallower than 
predicted from the accretion model. We interpret this 
as an indication of ram-pressure stripping; however, in 
Leo II, the cutoff is not nearly as drastic as Leo I, indi¬ 
cating that most of the chemical evolution had occurred 
before the MW began stripping and truncating Leo IPs 
SF. While this is by no means direct proof, the evidence 
supports the idea that Leo II chemically evolved without 
strong dynamical influence from the MW. 

Leo I on the other hand shows evidence of strong tidal 
interactions with the MW. The proper motions and dy¬ 
namic modeling of Leo I show it passed into the MW’s 
potential 2.33 Gyr ago, reaching i ts pericentric app roach 
1.05 Gyr ago (|Mateo et al.l 120011 iSohn et al.ll2ni3h . Ad¬ 


ditionally, the dynamic results correlate well with the 
CMD based SFHs. Leo I shows slow ancient SF, with an 
increase to its highest rate 2-3 Gyr ago, a fter which the 


galaxy stopped forming 

stars 1 Gvr asro (IGallaxt et al 

199S 

Gaouto et a 

.1119991 

Hernandez et al.||2000t Dolohin 

2002 

Weisz et al.1 

201411. 

Sohn et al.1 (1201311 suggest that 


the increased SF at 3 Gyr and the abrupt stop at 1 Gyr 
could have been caused by ram pressure compression or 
gravitational torques exerted by the MW as the galaxy 
passed into the MW’s potential. For Leo I, the dynami¬ 
cal history has had a significant impact on the SFH, and 
by extension on the chemical evolution. 

We interpret the sharp difference in the M parameters 
along with the dynamic information in the literature on 
Leo I, Leo H, IC 1613 and Phoenix, as possible evidence 
that the galaxies chemically evolved under different con¬ 
ditions. The first is passive chemical evolution, where 
gas rich, star forming galaxies (i.e. dir) gradually trun¬ 
cate their star formation as the gas supply is exhausted, 
evolving into a dTrans, producing a broad MDF in the 
process. This type of evolution is what is most likely 
what is seen in Leo H (dSph), Phoenix (dTrans) and 
IC 1613 (dir), at various stages in the process. Leo H 
formed the majority of its stars over 6 Gyr ago and the 
MDF suggests it chemically evolved without significant 
interactions before it was transformed into a dSph by the 
MW. Phoenix appears to be in the process of passively 
evolving; the dynamics and SFH show no indication of 
interaction. Instead the SFH shows a slow decrease in SF 
up to a few 100 Myr ago when SF stopped. Addition¬ 
ally, IC 1613 is isolated from both the MW and M31, 
and does not show dynamical evidence of interactions, 
and the MDF of IC 1613 is similar in shape to Leo H 
and Phoenix. It is known that dSph are created through 


the gravitational influence of a larger primary, if the in¬ 
teraction occurs concurrent with the chemical evolution 
it will be reflected in the MDF. Previous studies have 
shown that the MW’s influence on Leo I is reflected in 
the SFH and dynamics; we are suggesting it is also shown 
in the narrow MDF. 

7.2. MDF evolution 

One understated assumption in analytic chemical evo¬ 
lution models is that each model represents the end point 
of that galaxy’s evolution. However, dir and dTrans are 
not at the end of their chemical evolution since they 
are still forming stars and enriching the ISM. Therefore 
we must take care when comparing the MDFs of active 
galaxies such as IC 1613 and Phoenix, to the dSphs like 
Leo I and H, which arguably, have completed their chem¬ 
ical evolution 1 to 6 Gyr ago respectively. Additionally, 
for dSph the morphology-density relation suggest that 
the removal of gas interrupted the chemical evolution al¬ 
tering the chemistry to no longer reflect the end state of 
chemical evolution models. 

These various endpoints can be seen in our galaxy sam¬ 
ple. Leo I shows a MDF truncated on the metal rich side 
indicative of a halted chemical evolution. Leo H shows 
signs of mild truncation, as the GEM slightly overpre¬ 
dicts the number of metal rich stars and the accretion -|- 
truncation model was the best fit model for Leo H. As 
we discussed previously, Leo H is within the virial radius 
of the MW, but has a mostly circular orbit at a large dis¬ 
tance, therefore it has been weakly affected by the MWs 
gravity, as the mild truncation shows. 

As a dir, IC 1613 is still actively forming stars and en¬ 
riching its ISM, therefore we would not expect the MDF 
to reflect a fully completed chemical evolution. The on¬ 
going enrichment in IC 1613 manifests itself with fewer 
metal rich stars than predicted by the GEM. 

Phoenix, on the other hand, is the only galaxy to 
nearly match the metal rich end of the GEM. Phoenix 
is a transition galaxy, that expelled its last remaining 
gas reser ves with the final b urst of star formation ^ 100 
Myr ago (|Young et al.l 120071) . The isolated location and 
apparent exhaustion of gas imply that Phoenix has com¬ 
pleted its chemical evolution. 

One other consideration is that our field of view is small 
for all of the galaxies in our sample (covering 5-20% of the 
various galaxies). Given that many dwarf galaxies show 
metallicity gradients, and variations in SFH as a function 
of radius, the MDFs will also change as a function of 
radius, and the solutions to the OEMs cannot be applied 
to the entire galaxy. 

8 . CONCLUSION 

We measured the metallicity distribution functions of 
Leo I, Leo H, IC 1613 and Phoenix dwarf galaxies by 
measuring individual stellar metallicities, and by model¬ 
ing the stellar population to create synthetic MDEs for 
each galaxy. We find the synthetic MDFs to be a bet¬ 
ter representation of the overall metallicity distribution, 
because this method cuts down the photometric scatter 
propagated into the photometric metallicities. 

We ht each MDF with chemical evolution models. We 
find the simple model to be a poor fit for all of our dwarf 
galaxies, while the accretion -|- truncation model is the 
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best fit for the synthetic MDFs for all four of the galax¬ 
ies. The fact that all of our galaxies are best fit by the 
accretion -|- truncation model reflects the fact that most 
galaxies do not make it to the natural completion of the 
chemical evolution, either because they are still actively 
forming stars (dirs) or because they were truncated by 
other processes such as ram pressure stripping (dSphs). 

We find a similar accretion parameter, (M~4), for 
Leo II, IC 1613 and Phoenix despite the fact that they 
all have different masses, SFH, morphologies and aver¬ 
age metallicities. We interpret the resemblance in their 
MDFs as an indication that their chemical evolutions oc¬ 
curred under similar conditions, which indicates that Leo 
II completed most of its chemical evolution in isolation 
before it was significantly tidally disrupted and trans¬ 
formed into a dSph type galaxy by the MW. 

Leo I has the narrowest MDF, a much larger accretion 
parameter, and shows significant evidence for interac¬ 
tions with the MW at the same time as the galaxy was 
forming stars and chemically evolving. We suggest that 
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the MDFs can reveal dynamical interactions if they oc¬ 
cur in concert with SF and chemical evolution, or if the 
galaxy chemically evolved in relative isolation. The nar¬ 
rower MDFs are indicative of interactions shaping the 
galaxy’s present morphology, while a broader MDF indi¬ 
cates a passive evolution. The differences in the MDFs 
could be a way to distinguish between the two forma¬ 
tion pathways. To further this theory we would need to 
examine the MDFs and dynamic histories of many LG 
dwarf galaxies. 

Additionally, we measured metallicity gradients for 
Leo I and Leo 11. We see some evidence of gradients in 
Phoenix and IC 1613, however our data does not cover 
the radial extent required to determine a metallicity gra¬ 
dient greater than the error bars. 
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